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ABSTRACT 


The Geography Division of the U.S. Bureau of the Census must deter- 
mine the maximum potential extent of urbanized areas prior to conducting 
a Census of Population and Housing. Currently, the mapping task is 
cumbersome and time consuming, involving the manual integration of data 
obtained from many sources. As part of an ongoing research project for 
possible use in future censuses, new mapping approaches are being in- 
vestigated by the Jet Propulsion Laboratory. The Geography Division has 
noted that several of the map products may be quite useful in future 
mapping exercises. 

Anew technology for indicating potential expansion of urbanized 
areas has been developed. Through implementation of the system, Landsat 
imagery and more conventional types of cartographic data have been in- 
tegrated, and a data base has been produced from which various aspects 
of the urban environment can be represented. Thematic maps depicting 
areas, of urban expansion have been derived for three major urban regions 
in the United States: Orlando, Seattle, and Boston. Tabular reports 

summarizing those themes have been produced as well. 
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ft SECTION I 

t( 

THE CENSUS-URBANIZED AREA PROJECT 


A. INTRODUCTION , , ,, 

'• A o a V 

The primary focus .of this report is the research conducted at 
NASA/JPL in support of the Census-Urbanized Area Project (Reference 1-1), 
part of NASA's Application System Verification and Transfer (ASVT) pro- 
gram. Basically, JPL was given the responsibility of developing an 
information system capable of integrating data derived from Landsat 
imagery with data currently available to the Bureau of the Census. Also, 
JPL had been developing image processing techniques that could be uti- 
lized to map the urban fringe from Landsat Imagery. 

Most image processing for this project was completed by implement- 
ing programs from the Video Image Communication and Retrieval (VICAR) 
Image Processing System (Reference 1-2). VICAR was deyeloped in the 
early 1960s to support the imaging systems of JPL's unmanned planetary 
space program. Rather than a series of independent computer programs, 
VICAR is a fully-developed image processing system that includes a 
complete instruction set and a universal control card format. 

...... • . . 

The Image-Based Information System (IBIS) , a subset of the VICAR 
program library, was used extensively for data processing. Through 
implementation of IBIS, Landsat data and conventional data were inter- 
faced, and revised urbanized area maps were produced. In addition to 
revised urbanized area maps, statistical reports summarizing the results 
of the investigation were derived as well. 

o o 

B. BACKGROUND 

The U.S. Bureau of the Census has been mandated by the Federal 
Government to conduct a census of population and housing every five .. 
years beginning with the 1980 decennial census. As a consequence of 
this decision and associated legislation, the Census Bureau will have 
to process much more data during the next decade than in any previous 
decide. Furthermore, due to a steady increase in the population of the 
United States and a continued high mobility of the public, the enumera- 
tion process carried out by the Bureau of the Census will become more 
difficult. Accordingly, the Census Bureau is investigating means to 
revise existing practices to keep pace with the increasing demand for 
the collection and dissemination of data products. 

One of the most complex and more cumbersome of the procedures 
undertaken by the Geography Division of the U.S. Bureau of the Census is 
the mapping of the urban fringe prior to conducting each census. A 
variety of data must be collected and integrated to obtain maps of 
fringe areas. Conventional data collection arid integration practices 
have been found to be less than adequate for the derivation of urban 
fringe maps in a timely manner. Consequently, new technologies are being 
developed to avoid a potential backlog in obtaining map products. 
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In response to the needs of Che Geography Division, the Census- 
Urbanized Area Project was initiated. The purpose of this Project has 
been to test the utility of landsat data for the urban fringe mapping 
applications of the U.S. Bureau of the Census. Specifically, the pro- 
blems addressed in this ASVT are two: First, can Landsat data and image 

processing technology be used to obtain accurate renditions of urban 
fringe areas? Second, if such information is derived, can it be inter- 
faced with other materials available to the Census Bureau for/decision 
making purposes? Other factors such as the cost of data processing, the 
speed of information retrieval, the quality of output products, and the 
accuracy of those products are being reviewed as well. 


1. Data Processing Needs of the Bureau of the Census 

Beginning with the 1980 decennial census, the Geography Division 
of the Bureau of the Census must revise the urbanized area (UA) bound- 
aries of all Standard Metropolitan Statistical Areas (SMSA) in the U.S. 
every five years. While in 1970 there were 270 SMSAs, after the 1980 
census it is expected that there may be as many as 300 SMSAs. Since 
these. UA boundaries will be revised every five years, at least 600 UA 
maps will be constructed during the decade 1980-1989. To meet the needs 
rtf many Federal agencies that require knowledge of urbanized areas for 
decision making purposes, the Census Bureau will not be able to enjoy 
the luxury of processing UA maps on a weekly basis. If the mapping pro- 
ject is extended to cover the entire ten-year period, a minimum of 60 
maps must be updated each year. Based on that average, more than one 
map must be processed each week. 

a. Current Mapping Practices . Currently, urbanized area maps are 
compiled by integrating data obtained from a variety of sources. The 
updating process begins with a revised estimate of the extent of each 
urbanized area prior to conducting the census. The urbanized area bound- 
ary delimited as a result of the last census is integrated with informa- 
tion obtained from local and regional governments, aerial photographs, 
and population estimates and projections by the Census Bureau to define 
an updated UA boundary estimate. The region between the previous ur- 
banized area boundary ( inner line ) and the current estimate boundary 
( outer line ) is commonly referred to as the urban fringe zone (Figure 
1-1). Prior to conducting each census, the urban fringe zone is sub- 
divided into small geographic units or enumeration districts . After the 
census is conducted, statistics gathered for each enumeration district 
are analyzed, and decisions are made to determine if each enumeration 
district is urbanized. The decision is primarily based on population 
density. Finally, a new UA map is compiled by constructing a perimeter 
around the previous urbanized area and all newly-urbanized enumeration 
districts. 

The most time-consuming part of UA revision is the demarcation of 
the urban fringe zone. Although the. inner line is given, little is known 
about the position of the outer line. The process of combining the 
various data obtained from several sources is quite complex, and deci- 
sions are often subjectively based. Care must be taken not to omit any 
possible area from consideration. Consequently, the outer line of the 
urban fringe zone usually includes much rural land. If the outer line 
could be estimated more accurately, the subsequent census enumeration 





Figure 1-1. Dividing the Urban Fringe Into Enumeration Districts 

process could be streamlined .-.and, -as a result, the entire UA mapping 
project would become more ntf htgeabla. 

r 

b. Automated Mapping /j'P.yocedures . One purpose of the Census- 
Urbanized Area Project was to test the utility of Landsat imagery and 
image processing technology pi obtaining more accurate renditions of the 
urban fringe zone. Not only was it envisioned that the fringe area 
could be delimited more Accurately, it was believed that the entire 
process of/ estimating the urban fringe zone could be shortened with the 
implementation of computer technology. The use of an information system 
for interfacing remotely sensed data with other data available to the 
Geography Division of the Bureau of the Census was to be Investigated. 

In subsequent phases of the Project, remote sensing technology 
needed to derive urban fringe maps from Landsat data will be transferred 
to the. Geography Division. An information system used to integrate the 
remotely sensed data and the conventional data also will be provided. 

2. Principal Investigators 

Three principal investigating teams were selected to conduct re- 
search in support of the Census-Urbanized Area Project: (1) the Earth 

Resources Branch of the Goddard Space Flight Center (NASA/GSFC) (Refer- 
ences 1-3 and 1-4), (2) the Earth Resources Applications group at the 
Jet Propulsion Laboratory (NASA/JPL) (References 1-5, 1-6, and 1-7), and 
(3) the Image Processing and Analysis Center of the General Electric 
Company (GE) (References 1-8, 1-9, and 1-10). The primary responsibility 
of NASA/GSFC was coordination of the research efforts of the other two 
agencies. NASA/JPL was delegated the responsibility of developing a 
fully-automated information system capable of providing the products 
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needed to update urbanized area maps for the Geography Division. GE 
developed procedures required to transfer technology to the Bureau of 
the Census. GE also created a base-line mapping approach to the problem. 

The body of this report is divided into four sections. A descrip- 
tion of IBIS is covered in Section II, The urban fringe mapping tech- 
niques developed at JPL are covered in Section III. Description of 
three mapping applications are included in Section IV. Finally signif- 
icant conclusions drawn from the JPL effort will be summarized in 
Section V. A series of appendices that Include a technical description 
of IBIS, output products, and program timings are also included in this 
report. 


SECTION II 


THE IMAGE-BASED INFORMATION SYSTEM: 
AN OVERVIEW 


The Image-Based Information System (IBIS) Is a computer-based 
approach to the analysis of geographical situations* By grouping se- 
lected IBIS programs into processing operations, a variety of spatial 
phenomena can be investigated. The basic concepts of lata management 
and data processing within the Image-Based Information System are cov- 
ered in this section, | 

IBIS is considered to be a raster-baaed information system (Refer- 
ences 2-1, 2-2, 2-3). Most data entered into IBIS are in raster (image- 
based) format. However, the system is configured in such a manner that 
other data types, such as graphical and tabular data, may be used in 
analysis as well. 

Logical and mathematical interfaces have been provided to link 
all data files in an IBIS data-base superstructure (Figure 2-1). By 
utilizing these interfaces, information may be derived from simple 
associations of, or comparisons between, two or more data files stored 
in an IBIS data base. More complek procedures including polygon over- 
lay and cross-tabulation can also be investigated through IBIS techno- 
logy. 

A. DATA MANAGEMENT CONSIDERATIONS 

The image-formatted data plane is the primary data type utilized 
in IBIS processing. IBIS data planes may be obtained directly in image 
form, as in Landsat imagery, or they may be derived from data compiled 
by sources such as the U.S. Geological Survey, the U.S. Bureau of the 
Census, and the Defense Mapping Agency. „ 0 

Regardless of data type and origin, all data planes are incorpor- 
ated into a data set chat is referred to as the IBIS data base (Figure 
2-2). When investigating a specific problem, a data plane, may be in- 
cluded in, excluded from, or modified 1 'before any IBIS processing step. 

New data planes may be constructed with the system and may be used in 
subsequent processing steps. 


'The Image-Based Information System was designed by N. A, 
A. L r ,-^Zobris t of the Image Processing Laboratory at JPL. 
tion of IBIS included in this section is quite general, 
description of the system is- included in Appendix A. 
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Figure 2-1. A Configuration Diagram of the Image-Based Information System 
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Figure 2-2. Formation of an IBIS Data Base 

To maintain geometric consistency between all data planes includ- 
ed in the data base, an image plane exhibiting good radiometric or plani- 
metric qualities is designated the planimetric base . Digitized topo- 
graphic maps are excellent for this purpose, having been constructed 
from precise point datum and calibration. All other data planes are 
geometrically corrected to register with the planimetric base. 
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The us«?r of IBI^ can integrate various data types to forn&qn IBIS 
data base. Since the primary data structure is a raster format, image 
data planes are directly entered into the system; Graphical forms of 
data, Wmlly obtained in Cartesian reference form, must be transformed 
into im^ge space prior to inclusion as a data plane. Tabular data are 
not transformed into image space, but are linked to thd image data base 
through a logical interface. Data processing requirements for each data 
type are unique and will be covered individually. 

fl o ;j c 11 

a . Image Data . Most" image data sets entered into the data base 
are usually derived from Dandsat imagery or other multlspectral scanner 
sources. Other data are digitally encoded or scanned from aerial .photo- 
graphic products. Since image data planes are obtained from many 
sources, the spatial alignment of features contained in thosefrimages are 
often inconsistent from image to image. Consequently, provisions have 
been made to obtain a unified spatial surface through the implementation 
of resampling and spatial rectification procedures. 

Q 

Geometric correction procedures and spatial rectification routines 
for modifying image data are features of the VICAR image processing 
system. Although not specifically considered to be IBIS programs, these 
and other VICAR programs are necessary parts of the IBIS process. VICAR 
software is often used to obtain special information from Image data 
sources. Multlspectral classification of Landsat imagery is one example 
of a VICAR procedure othat: produces a data set used in IBIS procedures. 

b. G raphical Data . Graphical or vector data may also be entered 
into the IBIS data base. Graphical data are either produced locally on 
a coordinate digitizer, or are obtained from data tape. The Bureau of 
Census Urban Atlas file and DIME file are examples of data obtained on 
computer tapes.. Regardless of data origin, graphical data are trans- 
formed into image space prior to inclusion in the IBIS data base. 

As with image data, graphical data files must be in registry with 
the primary data base. Provisions have been made within IBIS to achieve 
the proper spatial alignment. (These corrections are made before the 
data are transformed into image space.) The registration technique" 
utilized is a two-step process. Initially, a general surface fit is 
achieved through the use of a least-squares affine transformation, then 
an exact geometric correspondence is obtained by implementing a spatial 
rectification procedure. The deformation of the original surface is 
controlled by the selection of tie points linking geographical features 
that are identifiable on both the graphical data file and the primary 
data base. 

Point and vector data are transformed into image space after geo- 
metric modifications are' completed. Three-dimensional or z-value data 
(x, y, z) are processed in a similar manner. The Cartesian reference 
components Of the data (x, y) are transformed into image Space coordi- 
nate values, while the z-value remains unchanged. 
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c . Tabular Data, Tabular data may be entered into IBIS via com- 
puter cards or digital tape. These data are stored in a tabular file 
that is linked to the data base through a logical interface. As a re-’J 
suit of this link, such data files are termed Interface files . 

2. The Georeference Plane 

One of the more important graphical data files entered onto the 
IBIS data base is the georeference plane . The georeference plane is a 
polygon file used in data aggregation and map generation procedures. 

In one case (the Census-Urbanized Area Project', the georeference plane 
has been constructed from census tract data obtained from the Urban 
Atlas files of the U.S. Bureau of the Census. 

Once the georeference plane is transformed into image space, each 
polygon, or region, must be identified. The region identification pro- 
cess involves the assignment of a unique value (or gray tone) to each 
individual region. After region identification, the georeference plane 
may be used in several higher-order IBIS procedures. For example, a 
polygon overlay of the georeference plane with some other image data 
plane can be Initiated (Figure 2-3), or the gr&y values of each polygon 
in the georeference base may be modified to produce a map depicting the 
results of a modeling application with data stored in an interface file. 

Several georeference planes may be included in an IBIS data base. 
For example, a data base may contain a census tract georeference plane 
and a congressional district georeference plane. The maximum number of 
regions that can be included in one georeference plane is virtually un- 
limited. Currently s up to 20,000 regions can be identified by gray 
value for an individual georeference plane. 

3. The Data Interface and Tabular Files 

All tabular files (interface files) are linked to at least one of 
the georeference planes included in the IBIS data base. The specific 
link is obtained by storing the numerical value (gray tone) representing 
each region of the georeference plane with tabular data describing at- 
tributes of that region (Figure 2-4). Attribute data may be statistical 
in origin, an identification code, or may be the. result of an image 
plane comparison routine such as polygon overlay or cross-tabulation. 

B. MANIPULATING DATA AND OUTPUT PRODUCTS 

If IBIS, or any other information system, were only a device to 
collect and store geographic data, the utility of the system would be 
quite limited. Most users of an information system require far more 
powerful features. Several means of data output, both pictorial and 
tabular, are needed. Also, the researcher may want to undertake complex 
modeling applications with the data files stored in the information 
system superstructure. 


CONCEPTUAL FRAMEWORK OF 
IBIS INTERFACE FILE GENERATION 
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Figure 2-3. The IBIS Polygon Overlay Procedure 
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Figure 2-4. The Date Interface 

Procedures for data output and data manipulation have been derived 
as part of the Image-Based Information System. Maps may be generated 
and tabular reports can be obtained. 

1. Data Manipulation Procedures 

Data stored in either the data base or an interface file can be 
modified or manipulated with IBIS software. New data planes and inter- 
face files are easily generated. Four basic data-manipulation proce- 
dures are currently available. 


a. Data Manipulation Betwee n Image Planes . New image data planes 
are generated as a function of two more image data planes. Chiefly, 
the procedures implemented to derive such data planes are VICAR routines, 
although some IBIS routines are also used. 

Simple transformations such as image addition, subtraction, multi- 
plication, and division are easily obtained. Complex functions are 
handled nearly as easily, and precise mathematical formulas may be 
specified. Image enhancement routines are available, as are several 
data classification and stratification routines. 

b. Data Manipulation Within the Interface File . Most functions 
available in the image domain are also available for analysis cl tabular 
data. Resultant from such operations, new tabular data entries are 
generated. Complex mathematical functions may be used to derive higher- 
order properties of data stored in an interface file. 

c. Data Manipulation of Image Data into Tabular Data . By imple- 
menting certain IBIS routines, data originally stored in image format 
may be summarized and copied into a tabular file. The majority of these 
routines are aggregation functions, an example of which is histogramming. 
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IBIS procedures for polygon overlay and cross-tabulation are within this 
realm of data transfer programs. 

d. Data Manipulation of Tabular Data into Image Data . The repre- 
sentation of tabular data in image form is primarily used as an output 
aid. By the implementation of a map generating routine, any georeference 
base can be modified as a function of an Interface file. Model ,ig of 
data is performed similarly. Data planes produced in this manner can be 
entered into the IBIS data base for subsequent operations. 


2. Data Output Features 

Two output formats are available to the system user: maps and 

tabular reports. Maps are produced directly from any image data plane 
or through modification of a georeference plane. Tabular reports are 
made available through the operation of a report generator. 


C. SUMMARY 

With a knowledge of image processing, an analyst can learn to 
operate the IBIS system. A researcher can utilize the system to store 
several data planes and much tabular data. With all of the information 
at the data user’s disposal, many complex modeling problems may be 
solved relatively effortlessly. 

The various modes of data entry, data manipulation, and data out- 
put provide the researcher with complete flexibility to structure a 
unique data base specifically designed for a particular problem or 
investigation. IBIS is merely a framework for analysis of spatial data 
The actual information system is constructed with the selection of 
specific image and tabular data. 
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SECTION III 

URBANIZED AREA MAPPING TECHNIQUES 


The Geography Division of the U.S, Bureau of the Census must fre- 
quently update urbanized area maps of all SMSAs in the United States. 

In response to a demand for more-current versions of these maps, the 
utility of new data sources and new technologies are being considered as 
a means to expedite the mapping process. One avenue of experimentation 
being investigated is the mapping of urbanized areas from Landsat imag- 
ery. 

A. A BASE FOR AUTOMATED MAPPING 

Landsat is an ideal (lata source for the production and revision of 
small-scale -maps depicting large urban regions. The areal coverage 
afforded by a single satellite image is quite extensive, and most urban 
regions are contained in a single frame. New imagery is made available 
on a regular basis. Consequently, map revision may be undertaken more 
frequently. Since imagery is available in series, maps depicting higher- 
order features of the urban environment, such as urban expansion, can 
be derived as well. 

1. Data Analysis 

To obtain urbanized area maps, a Landsat image must be analyzed 
and pertinent information must be extracted from the abundance of data 
contained in the satellite base. Analysis can be undertaken by visual 
techniques (much like visual interpretation of aerial photographs) , or 
with the aid of digital image processing. Visual analysis of Landsat 
imagery is an involved and time-consuming process that often requires 
the making of subjective decisions. For those reasons, visual analysis 
may at best provide an interim option for the Bureau of the Census. 
Conversely, digital image processing of Landsat imagery can provide 
similar information in a timely and objective manner. Once a specific 
technique has been developed, any SMSA can be mapped if satellite cover- 
age is available. Since a standardized sequence of processing steps la- 
being utilized, the analyst can be certain that the mapping of any SMSA 
will be a completely controlled process. 

2. A Methodology for Mapping the Urban Environment 

The Jet Propulsion Laboratory is investigating methods for the 
detection and mapping of urban land cover from Landsat multispectral 
scajiner (MSS) imagery. Through digital processing of Landsat data with 
the VICAR system, several land cover maps have been produced. Many of 
these maps are eventually used as image data planes in an IBIS data base, 
such as the IBIS modeling applications developed for this Project. 

Covered in this section is a description of a technique that has 
been developed to transform Landsat data into urban land cover maps. The 
mapping technique is referred to as a binary processor . A simplified 
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exp 1 of the binary processor is included in this section. More 

dot descriptions can be found in another report (Reference 3-1). 


3. The, Study Area 

As a means to further an understanding of the basic steps in the 
implementation of the binary processor, a mapping example has been in- 
cluded in this section: Landsat Imagery covering the Orlando SMSA has 

been selected, and land cover maps have been produced. The Bureau of 
the Census indicated that the Orlando region experienced a marked in- 
crease in population since 1970. It was hoped that the expansion of the 
urbanized area could be detected by this mapping technique. 

Two Landsat frames covering the Orlando SMSA were selected for 
analysis: (1) April 28, 1973 (1279-15285), and (2) April 18, 1975 

(1999-15091). Anniversary dates were selected to minimize any spectral 
variances that could be attributed to seasonal factors. 


B. PREPARATION OR LANDSAT DATA 

Several processing steps must be completed before urban land cover 
information is extracted from Landsat imagery. Initially, the imagery 
is obtained in the standard EROS data format (References 3-2 and 3-3). 

The imagery is reformatted to be compatible with VICAR system guidelines. 
The Landsat frame is modified to correspond to either a sinusoidal pro- 
jection or a "space oblique" Mercator projection (Reference 3-4) . 

Following reformatting, the region containing the study area is 
extracted (Figure 3-1). At the analyst's discretion, cosmetic procedures 
such as dcstriping and bad-line removal are performed. 


1. Planimetric Modifications 

Several Landsat images are processed to obtain maps depicting 
changes in land cover or other dynamic aspects of the urban environment. 
To facilitate comparisons between maps or images, a planimetric data base 
is selected for the purpose of registering these data planes. A plani- 
metric base may be a topographic map, a UTM grid, or any other geographic 
expression of the Earth's surface. Often, a Landsat image can be selected 
as the planimetric base when the study area is of a limited size, and 
precise radiometric accuracies are not critical, 

The Landsat images are registered with the planimetric base 
through the implementation of a surface-fitting algorithm. Tie points 
selected to define the deformational characteristics of the surface are 
found with the aid of an interactive duoimage pattern recognition and 
comparison routine. An equation of transformation is derived, and the 
image is spatially rectified to coincide with the planimetric base. 





Figure 3-1. A Landsat Image Covering the Orlando SMSA 


lnformat ion Enhancements 


After all needed imagery is registered to the planlmetric base, 
the next phase of data preparation is begun. The raw satellite data is 
transformed into a more usable data plane (Figure 3-2) through the imple- 
mentation of an information-enhancement routine. Simply, the process of 
information enhancement yields an image that contains more usable infor- 
mation for the analyst. Ratio (Reference 3-. t >) and eigenvector (principal 
components) (Reference 3-6) transformations are routines commonly used 
for information enhancement. 


In the case of the Census-Urbanized Aren Project, information- 
enhancement routines were implemented to enable easier distinction be- 
tween urban and nonurban land cover- In another research project, it 
was found that ratio-transformed images were especially useful in the 
mapping of land cover from Landsat imagery for areas of Los Angeles, 
California, and Madison, Wisconsin (Reference 3-1). A similar technique 
(Reference 3-7) was used to investigate land-use change around Atlanta, 
Georgia. Later investigations into eigenvector processing have also 
proven quite promising for improving data quality. 

One consequence of the utilization of information-enhancement 
routines was the reduction in the number of data planes needed for 
mapping the urban landscape. When as many as four raw-data planes were 
required for classification of land cover with conventional operations, 
only one data plane was required to achieve the same classification 
after an information-enhancement routine had been implemented. With a 
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Figure 3-2. The Second Data Vector of an Eigenvector 
Transformation Implemented to Enhance Information Content 


reduction in the number of data planes needed tor classification, a 
complex Bayesian or maximum- likelihood classifier was no longer needed. 

A simpler and much faster classifier could he utilized to derive similar 
results . 

C. THE BINARY PROCESSOR 

The Geography Division of the U.S. Bureau of the Census is pri- 
marily concerned with distinguishing urban from nonurban land cover. To 
obtain such information, the classification of Landsat imagery was based 
on a binary decision: Is the land urban or nonurban? A simple binary 

decision processor was implemented to obtain the desired maps. 

1. Theoretical Considerations 

Binary classification of land cover is based on analysis of a 
histogram sampled from a Landsat image covering a study area. Simply, 
a histogram is considered to be a spectral s i gnat ure characterizing the 
spectral attributes of all features contained in the satellite image. 
Given that the implementation of an information processing algorithm has 
actually facilitated the distinction between urban and nonurban land 
cover, a histogram sampled from that image should also indicate the dis- 
tinction. More specifically, the histogram consists of two predominant 
spectral signatures representing urban and nonurban land cover. Further- 
more, the histogram should tend to be a bimodal distribution (Figure 3-3) 
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Figure 3-3. Two Spectral Signatures Derived From 
a Bimodal Distribution Sampled From a Landsat Image 

a. Threshold Analysis . Through closer inspection of the histo- 
gram, an intersection between the two spectraL signatures can be found. 
The point of intersection is referred to as the t hreshold limit of the 
distribution. The threshold limit separates the distribution into two 
distinct ranges of data. Each range corresponds to the unique spectral 
signature of urban or nonurban land cover, 

b. Classification . After determining the threshold limit, the 
Landsat image that has been the source of data sampling can be classi- 
fied. As with the histogram, the data contained in the Landsat image 
are divided into two ranges. Data in one range correspond to urban 
land cover. When each range of data is symbolized, a two-class land 
cover map (Figure 3-4) is produced. 

2. Application of the Binary Processor 

Two mapping exercises were completed by processing imagery 
covering the Orlando area. Each exercise was based on the implementa- 
tion of the binary processor. The primary difference between the two 
exercises was that a different information enhancement routine was 
utilized in the course of completing each map. 

a. An Application Based on a Ratio Transformation. Application 1 
Previous research indicated that effective land cover maps could be 
derived from ratio-transformed data. Consequently, such imagery was 
utilized in the first mapping application. Usually, the ratio trans- 
formation used in analysis is formed when the red data band (MSS 5) is 
divided by the first infrared data band (MSS 6) of the Landsat MSS. 
However, in the Orlando case study, the data were useless due to sevefe 
striping of MSS band 6. Instead, MSS band 7, the other infrared data, 
band, was utilized in the ratio transformation. 
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Figure 3-4. Diagrammatic Depiction of the Binary Processor 

A histogram (Figure 3-5) was generated by sampling a portion of 
the transformed image. Through visual analysis of the histogram, two 
distinct spectral signatures representing nonurban and urban land cover 
were identified. A threshold limit was determined, dividing the distri- 
bution into two data ranges. In previous research (Reference 3-1) it 
had been determined that the urban spectral signature (characterized by 
light materials such as concrete) was confined to the upper spectral 
range, and nonurban to the lower range (where dark-colored materials 
were predominant). 

Based on the analysis of the histogram, the ratio-transformed image 
was divided into two ranges of data. Each range was symbolized to de- 
rive an ur' ,i. land cover distribution map (Figure 3-6). 

FREQUENCY DISTRIBUTION 



Figure 3-5. A Histogram Used in the Analysis of the 
Ratio Transformation: MSS Bands 5/7 


3-6 





Figure 3-6. The Binary Classification of the Ratio 
Transformation. Features Depicted in Black are Urban 

The distinction between urban and nonurban land cover was not 
particularly effective^. Although major urban regions were identified, 
a significant amount of nonurban land was also erroneously included in 
the urban class. Consequently, classification based on the ratio trans- 
formation was not useful in this application. 

b . An Application Based on Eigenvector Data, Application 2 . 

A second attempt to map the distribution of urban land cover around 
Orlando was initiated. In the second application, the raw data were 
transformed by the implementation of an eigenvector algorithm. Of the 
four data vectors derived, data vectors 1 and 2 (F'p.ure 3-7) were used 
in the classification of land cover. 

In the previous mapping exercise, one iteration of the binary 
processor was utilized to obtain a land cover map. However, two itera- 
tions were required when eigenvector data were utilized. One iteration 
of the binary processor was implemented to remov iny conflicts between 
the spectral signatures of urban land and water. 


A major cause of the misclassif ication has been attributed to signal 
attenuation caused by water satuiated soils. The signal attenuated 
tends to modify the nonurban spectral signature in such a manner as 
to reseaibie urban spectral signatures. 
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Figure 3-7. The Second Data Vector of an Eigenvector 
1 ransformation 


Histograms (Figure 3-8) were sampled from eigenvectors 1 and 2. 
The distribution of eigenvector 1 was used to separate land from water 
and a map depicting those features was produced. Hie water features 
mapped from vector 1 were subtracted from the second eigenvector image 
removing any conflicts between the spectral signatures of water and ur 
ban land. Then a histogram was sampled from the second data vector, 
and a map depicting the distribution of urban land cover was produced 
in the usual manner. 


Since the Geography Division was interested in sequential mapping 
for updating purposes, two Landsat images (1973 and 1975) were classi- 
fied (Figures 3-9 and 3-10). Because the raw-data images were regis- 
tered to the planimetric base, all resultant mans were registered to the 
base as well. Consequent ly , maps depicting, changes in the distribution 
of urban land cover could be easily obtained. Urban expansion between 
1973 and 1975 has been documented in this manner (Figure 3-11). 


SUMMARY 


The procedures utilized for the preparation and classification of 
Landsat imagery to obtain land cover maps have been covered in this 
section. As a result of the analysis, a map depicting the distribution 
of urban land in 1975 (see Figure 3-10) has been produced. That image 
has been utilized as an image data plane in the IBIS mapping application 
covered in the next section of the report. 
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Figure 3-8. Two Histograms Used in the Analysis of Eigenvector 


Figure 3-9. The Binary Classification of the Second Data Vector of 
the Eigenvector Transformation, 1973 





Figure 3-10. The Binary Classification of the Second Data Vector of 
the Eigenvector Transformation, 1975 
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SECTION IV 


TECHNOLOGY APPLICATIONS 


The Jet Propulsion Laboratory has developed a complete mapping 
system designed to meet the needs of the Geography Division of the U.S. 
Bureau of the Census. The mapping system has been designed for the pur- 
pose of updating urbanized area boundaries on maps covering major urban 
regions in the United States. To achieve that purpose, procedures need- 
ed for data-base management, urban land cover mapping, and thematic map 
generation have been designed as basic components in the system. 

The mapping system is based entirely on routines contained in 
either the VICAR image processing system or the Image-Based Information 
System (IBIS). The basic features of both VICAR and IBIS used in the 
mapping system have been described in previous sections of this report. 
In this section, application tests using the mapping system wilJ;/be 
described. // 

0 

Three applications of the VICAR/IBIS-based mapping system were 
undertaken. These applications were specifically devised for the pur- 
pose of testing the system design and verifying the utility of the 
system for mapping urbanized areas. To fully demonstrate the complete 
features of the JPL approach, the utilization of different processing 
schemes was emphasized in each application. 

In the first mapping exercise, urbanized areas within the Orlando 
SMSA were mapped. As the initial test of the system, basic data-base 
management procedures and mapping techniques were developed. The next 
test was implemented on data covering a portion of the Seattle-Everett 
SMSA. Emphasis was placed on the ability to utilize land cover data 
compiled by other researchers. The final mapping application covering 
portions of seven SMSAs around Boston, Massachusetts, was the most 
complex test of the mapping system. In that application, data planes 
were completed only when several smaller data planes were combined. 

A. MAPPING URBAN EXPANSION WITHIN THE ORLANDO SMSA 

The population of the Orlando SMSA has substantially increased 
since the 1970 census. It is highly probable that urbanized regions 
within the SMSA have also expanded during that period, and, as a result, 
the urbanized area map compiled from 1970 census data required revision. 
Through the implementation of VICAR and IBIS programs, information per- 
taining to urban expansion was gathered for map revision. 

1. Preparing the Data Base 

Derivation of the Orlando urban expansion map required the crea- 
tion of a data base. The data planes comprising the data base were ob- 
tained from many different source'-, since complete information was not 
available from any one source. Two primary data planes directly entered 
into the data base were a census tract outline and a Landsat image. 
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a. Selection of the Planlmetric Base 
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Figure 4-1. The Planimetric Base for the Orlando SMSA 
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data planes Including another Landaat image, a census tract outline 
image, and all other data planes derived through IBIS processing, were 
registered to the planlmetric base. 

b. Forming the Georeference Plane . Another important data plane 
in the data base was the georeference plane. Because of the Geography 
Division's familiarity with census tract data, a logical choice for the 
georeference image was a census tract map. The georeference plane, 
based on census outline data, was not directly entered into the data 
base. Rather, it had to be generated through IBIS processing steps. 

4 

Initially, the census tract data were obtained in Cartesian refer- 
ence form and had to be translated into image-based coordinates. Then 
the coordinate data were modified to correspond to the surface character- 
istics of the planlmetric base. First, a least-squares surface fitting 
algorithm was implemented to obtain a general fit, after which a spatial 
rectification algorithm was invoked to remove any local distortions. 

Once the desired surface correspondence was achieved, the image-based 
coordinates were transformed into image space, resulting in the genera- 
tion of another data plane (Figure 4-2) . To verify the accuracy of the 
surface fit, the census tract image and the planlmetric base could be 
combined (Figure 4-3) . 

Once the census tract data was transformed into image space, the 
actual georeference plane was created by assigning each census tract a 
unique value (gray tone). In the Orlando SMSA, 101 census tracts were 
identified in this manner (Figure 4-4). It should be noted that because 
the georeference plane had been derived directly from the census tract 
outline data plane, the georeference plane was already registered to the 
planimetric base. 

To establish a link between a georeference plane and other forms 
of data, an interface file was produced. The interface file is a tabu- 
lar type of data file designed to contain pertinent information about 
each region in the georeference plane. Besides containing the linking 
information (gray value codes) for each census tract, the interface file 
commonly contains other types of data such as population counts, area 
calculations, and land cover aggregations. 

2. Preparing the Thematic Data Planes 

In the Orlando research application, a decision was made to map 
urban expansion for the period 1970-1975. Two data planes were inte- 
grated to produce a map depicting that theme. One data plane, derived 


Census tract data used in the Orlando application was digitized by the 
Geography Division of the U.S. Bureau of the Census. The data, stored 
with data for other SMSAs, was accessed from a digital tape referred 
to as the Urban Atlas File tape. The data files, including the Orlando 
data file, contain a record of census tract identification codes and 
digitized coordinates outlining each census tract in the SMSA. 
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Figure 4-2. The Census Tract Data Plane for the Orlando SMSA 

from census data, represented urbanized areas in 1970. Another data 
plane, obtained through analysis of Landsat imagery, portrayed urbanized 
lands for the later date. 

a. Forming the 1970 Urbanized Area Map . Population density is 
the primary factor considered by the Geography Division when determining 
if an area is urbanized. Generally, if an area has a population density 
of at least 1000 people per square mile, the area is considered to be 
urbanized. While other factors are also considered, they are not of 
major concern for this specific application. A map of urbanized areas 
is produced by delineating all areas that have been designated as urban, 
based on that criterion. Usually, the urbanized area map is established 
from data gathered at the block level. 

At the time the Orlando research project was undertaken, the 1970 
urbanized area map produced by the Bureau of the Census was not available. 
Consequently, an urbanized area map (Figure 4-5) was computed through 
IBIS processing. Simply, population density values were computed with 
population and areal data previously stored in an interface file. Then, 
based on the 1,000-people-per-square-mile criterion for urbanized, land, 
two classes of data were derived. Finally, a map-generating routine 
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Figure 4-3. The Integration of the Two Data Planes to Check 
Registration Accuracy: the Orlando SMSA 


was implemented to modify gray va lues of each region in the census tract 
georeference plane by a function of that data division. The resulting 
map depicted urbanized regions for the Orlando SMSA in 1970. That 
urbanized area map, established from data gathered at the census tract 
level, was not finely detailed when compared to the typical urbanized 
area maps compiled by the Bureau of the Census. However, the map did 
present a useful reference point for delineating urban growth. 

b. F orming the 1975 Urbanized Area Map . A map typifying the 
distribution of urbanized land in the Orlando SMSA during 1975 was 
derived from Landsat imagery (1999-15091, May 1975). The Landsat 
imagery were also data planes in the Orlando data base. Before an urban- 
ized area map (Figure 4-6) could be generated, the imagery had to be 
transformed to increase information content, and then a binary processor 
was implemented to classify the data. Those procedures required for 
information enhancement and classification have been detailed in Section 
3. 
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Figure 4-4. The Georeference Plane for the Orlando SMSA 
3. Synthesis of the Data Planes 

Once the two thematic data planes representing urbanized lands in 
1970 and 1975 were generated, a new thematic data plane (Figure 4-7) 
portraying urban expansion between 1970 and 1975 could be produced. The 
new thematic data plane was derived through the synthesis of the two 
initial urbanized data planes. 


a. V isual Acuity . One consequence of integrating data obtained 
from two different sources into a single map form is that differing 
spatial characteristics of the data become evident. The 1970 data, 
being a result ot areal generalization, were delineated along census 
tract boundaries and appe.*-ed as a highly simplified feature. Tine 1975 
data, obtained from Landsat imagery, were of a much finer resolution and 
appeared as a highly complex feature. 


The nature of data representation is chiefly a function of differ- 
ing modes of data collection and symbolization. Population statistics 
and areal measurements used in the derivation of the 1970 urbanized area 
map were accumulated by census tract. Because calculated population 
density values used to define the urban/nonurban nature of a tract do 
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Figure 4-5. The 1970 Urbanized Area Map of the Orlando SMSA 
Derived From Population Density Statistics 

not indicate where the urbanized areas are distributed within a tract, 
the only meaningful representation of urban land in 1970 was at the 
tract level. Conversely, the map depicting urbanized areas in 1975 was 
based on the analysis of Landsat imagery. Landsat data are collected at 
the much finer resolution of an 80 x 60 meter grid. The daca are not 
aggregated before analysis, and can be depicted at that fine resolution 
in the output product. 

b. Statistical Reports . In addition to the 1970-1975 urban ex- 
pansion map, a statistical report^ of urban expansion was obtained as 
well. The report was derived through the implementation of various IBIS 
procedures. First, the land cover counts for each type of land cover 
(nonurban, urban in 1970, and urban by 1975) were tabulated by census 
tract. The tabulation was obtained through polygon overlay of the the- 
matic map and the georeference plane (of the census tracts), then per- 
centage coverages for each type of land cover were computed by census 
tract, and finally a report was produced via a report generator. 


All statistical reports are included in Appendix B. 



The 1975 Urbanized Area Map of the Orlando SMSA Derive 
from Landsat Imagery 
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Urban Expansion in the Orlando SMSA for the Period 1970 
1975 




1. 


Preparing the Data Base 


a. Selection of the Planimetric Base. Because a thematic 


land cover classification of Landsat imagery was available, the Landsat 
image from which the classification was created was a logical choice for 
the planimetric base (Figure 4-8). A Landsat image (1690-18245, June 
1974) and the classified scene derived from a 1975 Landsat image (2 182— 
18201, July 1975) were precisely registered to a map base where true 
north was directly vertical. Rotation due to satellite trajectory was 
removed and local surface distortions were minimized. As in VICAR and 
IBIS processing, the imagery was precisely registered to a map base with 
the implementation of a surface fitting routine. 


4-9 





ORIGINAL PAGf IS 
OF POOR 


4-11 


PRECEDING PAGE CLANK NOT FILMED 




b. Forming the Georeference Plane . The georeference plane was 
prepared in the usual manner. Census tract data, obtained from an Urban 
Atlas file, were registered to the planlmetrlc base, and a census tract 
data plane** (Figure 4-9) was created. The census tracts were then 
assigned unique values (gray tones) tv form the georeference plane. 
Finally, an interface file was generated. 



Figure 4-9. The Census Tract Data Plane Overlayed on the Planimetric 
Base of the Seatt le-Everett SMSA 


It should be noted that sone census tracts in the eastern portion of the 
SMSA were deleted from the georeference base. They were not included 
because satellite imagery is not available for those areas. 
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2 . 


Preparing the Thematic Data Planes 


a. Forming the 1970 Urbanized Area Map . As in the Orlando 
experiment, the 1970 urbanized area map (Figure 4-10) was derived thtough 
a calculation of population density. The population statistics, obtained 
from the 1970 census, were divided by IBIS-calculated area measurements 
for each census tract. 
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Figure 4-10. The 1970 Urbanized Area Map of the Seattle-Everett SMSA 
Derived From Population Density Statistics 


b. Forming the 1975 Urbanized Area Map . Instead of a JPL 
classification of the land cover, the land cover classification of a 
Landsat image recently completed by the USGS was utilized. Originally, 
twenty-one classes were included in the map. However, that level of 
detail was not needed for the Seattle-Everett application. The number of 
classes were reduced to three (urban, nonurban, and water) before the 
data plane was finalized (Figure 4-11). The pattern of urban land cover, 
the darkest tone, is clearly visible. 
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Figure 4-11. The 1975 Urbanized Area Map of the Seattle-Everett SMSA 
Derived From Landsat Imagery 

3. Synthesis of the Data Planes 

Once the 1970 and 1975 urbanized area maps were completed, they 
were combined as in the previous experiment. The result of the synthesis 
was a map depicting urban expansion between 1970 and 1975 (Figure 4-12). 

The visual appearance of the map is quite similar to that for 
the Orlando SMSA. The 1970 distribution, based on areal measurements, 
is very generalized. The 1975 urbanized area additions are point-specific 
data derived from the more-finely resolvable Landsat data source. 

The USGS is still in a research mode for the generation of these 
andsat-based land cover maps. However, when they are made available, 
the utilization of these maps as a data plane will greatly reduce the 
processing time required to generate urban expansion maps. Other data 
collecting and mapping agencies may also be able to supplement the 
Bureau of the Census with land cover maps for major urban regions. 
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Figure 4-12. Urban Expansion in the Seatt le-Everett SMSA 
for the Period 1970-1975 
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MAPPING URBAN EXPANSION WITHIN THE BOSTON REGION 


For urban area analysis, the Bureau of the Census often requires 
precise delineation of the urban fringe. Since the utilization of popu- 
lation density values measured by census tract can only be at best a 
rough approximation of the urbanized ^trea boundary, the Geography 
Division must consider other, more detailed, geographic parameters for 
determining the urban fringe. For the application covering the Boston 
region, the actual 1970 urbanized area boundary was digitized and pro- 
cessed. The boston application represented the most complex of situa- 
tions, as seven SMSAs were included in the study region. 

Instead of providing maps of urban expansion for final products 
as in the previous experiments, several new products were produced for 
the Boston application. The products provided the Bureau of the Census 
with a variety of ancillary photoproducts with increased spatial acuity. 
The products might be included with the digitally-processed urban expan- 
sion data to enhance the locational factors of those maps. The products 
could also be used within the more traditional framework of photointer- 
pretation practiced by the Geography Division to provide a Regional base 
for their mapping applications. 


Preparing the Data Base 


a. Selection of the Planimetric Base . As in previous experiments, 
Landsat imagery was selected for the planimetric base. However, two 
Landsat scenes (20829-14311 and 20829-14317, both imaged on April 30, 

1977) were combined to form this planimetric base. The Boston study 
region, covering an extremely large area, was contained in two sequential 
images. The images were intersected to obtain the desired study region 
(Figure 4-13). 


! i 


In an attempt to assess the effect of new data formats with the 
upcoming Landsat 3 im^ery, the scale of resolution for the planimetric 
base was modified. Instead of the normal 60 x 80 meter resolution, the 
imagery was resampled to an effective resolution of 57 x 57 meters. No 
major degradations were initially apparent as a result of the resampling. 

b. Forming the Georeference Plane . Portions of the seven SMSA re- 
gions were included in the Boston study, region. These SMSAs were indi- 
vidually processed to form the complete georeference plane. Each SMSA 
was registered to the planimetric base individually, since distortions 
were a localized problem. Six complete SMSAs were included in the data 
base: (1) Boston, (2) Brocton, (3) Fitchburg-Leominster , (4) Lawrence- 

Haverhill, (5) Lowell, and (6) Worcester. A portion of the Providence 
SMSA was also included. 




x 
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At the time of the 1970 census, some portions of the Boston region 
were not included in the network of SMSAs along the eastern seaboard. 
Census tracts and municipal boundaries were added in those areas to 
complete the georeference plane. This ancillary data was digitized from 
maps prepared by the Geography Division. Clearly, the census tract data 
plane (Figure 4-14) for the Boston area is very complex. .'-X 
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Figure 4-14. The Census Tract Data Plane for the 
Boston Region 

c. Forming the 1970 Urbanized Area Map . Instead of utilizing 
population density values to determine urbanized areas within the Boston 
region in 1970, a new data plane was Included in the data base. The 1970 
urbanized area data plane was based on the inner line map developed by 
the Bureau of the Census after the 1970 census. All regions within the 
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Inner line were considered to be urbanized In 1970. Through analysis of 
the 1970 urbanized area data plane (Figure 4-15) It was noted that ur- 
banized area boundaries did not always conform to specific census tract 
boundaries . 



Figure 4-15. The 1970 Urbanized Area Map 
of the Boston Region 


The final product from the Boston application was intended for 
manual photointerpretation of urban expansion between 1970-1977. The 
1977 Landsat image used to form the planimetric base was modified to 
highlight the areas urbanized by 1970 with an orange tint (Figure 4-16). 
With a typical color reconstruction of a Landsat image, urban regions 
are characteristically represented by d blue tones. With the 
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Figure 4-16. A Color Reconstruction of a Landsat Image 
of the Boston Region Where the Urbanized 
Area in 1970 has been Modified 
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integration of the 1970 urbanized area information* only those urban 
areas that were not considered urbanized in 1970 retained a true color 
representation, allowing immediate recognition of urban expansion since 
1970. ^ 

D. SUMMARY 

The research covered in this section represents the culmination 
of JPL's research in support of the Census-Urbanized Area Project. 

Through the utilization of the mapping system developed at JPL, the needs 
of the Geography Division of the Census Bureau can be satisfied. 

Landsat imagery may be quickly classified into urban and nonurban 
classes. Data for 1970 can be derived as well. The various data planes 
can be combined to fora the map depicting urban area expansion, and 
statistical reports summarizing urban land cover changes can also be ob- 
tained. These maps may be displayed at any desired map scale, because 
the photorecorded image is easily enlarged. Consequently, the Bureau 
of the Census geographer can assess the need for revision of enumeration 
district boundaries for the upcoming census. 
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SECTION V 


CONCLUSIONS 


A mapping system has ‘been developed by the Jet Propulsion Labora- 
tory to serve the needs of the Geography Division of the U.S. Bureau of 
the Census. Through implementation of VICAR and IBIS programs, the 
basic components of the system, the Geography Division can obtain maps 
depicting many characteristics of the urban environment. By integrating 
Landsat imagery with more traditional data types, elements of urban 
expansion are easily identified. The information depicted on such maps 
could be useful for identification of urbanized area boundaries and re- 
vision or addition of enumeration districts prior to the 1985 census. 

A. SYSTEM FLEXIBILITY 

Depending on guidelines formulated by the Geography Division for 
the Census-Urbanized Area Project,' a sequence of image processing steps 
can be designed to produce specific maps or products requested. A 
variety of such products have been presented in this report. 
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1. Levels of Technology 


In the body of the report, many maps and map-like products have 
been utilized to represent intermediate steps In mapping processes. 
Consequently, the attributes of these maps have not been emphasized. 
However, in some instances there may be sufficient mapped information 
to be useful in updating urbanized area boundaries. IBIS processes can 
be categorized within a hierarchical structure based on technology for- 
mulated for a previous level or levels. While each higher level is 
based on more sophisticated data processing procedures, the time required 
for map analysis after data) processing is shortened with each techno- 
logical advancement. 


a. Level 1: i ;e Enhancement . The most basic approach to the 

problem of revising urban ! ' area maps is through visual analysis of 

color-enhanced Landsat imag As exemplified in the Boston applica- 

tion, sufficient information can be obtained from such imagery (see 
Figure 4-13) to provide a trained geographer with the data needed to 
identify urban regions. When the 1970 urbanized area mask is integrated 
with the satellite scene (see Figure 4-16), assessments of urban expan- 
sion can also be made. 

b. Level 2: Addition of a Georeference Plane . The addition of 

a georeference plane such as a census tract boundary image (see Figures 
4-9 and 4-15) can simplify the task of urban area analysis. Since the 
visual stimulus of the census tract data is easily assimilated by the 
geographer, the spatial aspects of the satellite-image phenomena can be 
more clearly identified. 

c. Level 3: Thematic Maps . A more sophisticated approach to 

urbanized area mapping is based on the utilization of land cover maps 
(see Figure 4-11) produced through multispectral classification of 
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POPULATION STATISTICS OiTADCB r»» THt 1*70 Cf'lSUS OF POPULrtTHU AND HOUSING 

Figure 5-1. Choroplethic Map of Population Density of the Orlando SMSA 


Land sat imagery. Tills analytic approach to data collection provides an 
objective basis for future map revision applications. 


d. Level 4: Data Int egration. The most elaborate means to re- 

vise urbanized area maps is based on the integration oi conventional 
data types (such as population statistics) with Landsat imagery. Infor- 
mation that can not be derived from individual data sources can be ob- 
tained after the two data planes are integrated. Urban expansion be- 
tween 1970 and 1975 has been mapped by this process for the Orlando 
SMSA (see Figure 4-7) and the Seattle-Kverelt SMSA (see Figure 4-12). 


2. Other Thematic Maps 

A statistical map generation program has been developed for the 
construction of thematic maps from data stored in IBIS data files. 

Through implementation of this program, a choroplethic map (Figure 5-1) 
depicting levels of population density in 1970 has been produced for the 
Orlando SMSA. All data required for the production if this map were 
included in the interface file and georeference plane used on the Orlando 
appl icat ion. 
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The mapping system Is easily adaptable to atlas generation. 
Statistical distributions can be calculated with data stored in an inter- 
face file, and the georeference plane can be modified to depict the data 
distribution. Gray tones and patterns are easily selected for symboli- 
zation of spatial data. 

B. COST OF OPERATION 

1. Main-Frame System 

The Geography Division is naturally concerned with the expected 
costs of operating any mapping system where Landsat data is analyzed in 
a digital mode. Selected program timings for the three mapping applica- 
tions completed during phase I of the Program are included in Appendix 

C. It is made evident in the Orlando case, where IBIS and non-IBIS pro- 
cessing times are Itemized, that IBIS data processing does not signifi- 
cantly Increase the computer time expended in the mapping process. The 
added benefits of the products derived through IBIS software should 
outweigh any increased costs. 

2. Minicomputer Systems 

All data processing described was performed on an IBM 360-65 
computer located in the Image Processing Laboratory (IPL) at the Jet 
Propulsion Laboratory. Since the Geography Division may have difficulty 
in accessing a large, main-frame computer under normal operating condi- 
tions, the options that utilize a minicomputer system may be preferable. 
During 1979, JPL completed mimicomputer versions of VICAR and IBIS. 

These versions of "mini-VICAR" and "mini-IBIS" can be made available to 
the Bureau of the Census or other potential system users upon a request 
made to NASA-COSMIC. 

Although testing of the mapping system with the minicomputer 
system has not begun, three limiting factors about system characteristics 
are known. First, since minicomputers do not have large direct-access 
memories, many data sets will reside on peripheral storage units. Con- 
sequently, input/output timings will be increased as these files are 
transferred to and from memory. Second, certain programs that are 
' compute^bound ' will require more time for completion because the central 
processing units of most minicomputers are slowe/T than main-frame de- 
vices. The spatial rectification and region-identification routines will 
be most affected. Third, there will be certain limitations to data set 
sizes. Some large regions may need to be segmented for processing in a 
minicomputer environment. 

The minicomputer approach offers several processing advantages. 

The hands-on capabilities of minicomputer systems will enable faster 
processing of data. Routines that require several iterations to achieve 
desired results will be completed more quickly, as errors are identified 
and corrections are made during data processing sessions. All inter- 
mediate steps in a mapping procedure can be monitored. Consequently, 
back tracking usually caused by the unavailability of intermediate pro- 
ducts can be eliminated in the minicomputer mode. 


5-3 


REFERENCES 


1-1. NASA, Project Plan: Applications Systems Verification Test (ASVT) . 

Census Urbanized Area , NASA Goddard Space Flight Center, Greenbelt, 
Maryland, 1979. 

/( 

1-2. Seidman, J. B., and A. Y. Smith, VICAR Image Processing System : 

Guide to System Use . JPL Publication 77-37 Revision 1, Jet Propul- 
sion Laboratory, Pasadena, California, 1979. 

1-3. Christenson, J. W., J. B. Davis, V. J. Gregg, H. M. Lachowski, and 
R. L. McKinney, Landsat Urban Area Delineation , Intralab Project 
75-3, NASA/Goddard Space Flight Center, Greenbelt, Maryland, 1977. 

1-4. Christensbn, J. W. , et al . , "Urbanized Area Analysis Using Landsat 
Data," paper presented at the conference on Computer Mapping, 

Harvard University, Cambridge, Mass. , 1978. 

1-5. Davis, B. , and S. Z. Friedman, "Assessing Urbanized Area Expansion 
Through the Integration of Landsat and Conventional Data," 

Proceedings, 45th Meeting of the American Society of Photogrammetry , 
pp. 776-791, 1979. 

1-6. Friedman, S. Z., and G. L. Angelici, "The Detection of Urban Ex- 
pansion from Landsat Imagery," Remote Sensing Quarterly , Vol, 1, 

No. 1, 1979; presented at Annual Spring Meeting of the Association 
of American Geographers, 1978. 

1-7. Wherry, D. B., and S. Z. Friedman, "Cartographic Applications of 
an Image Based Information System," Proceedings, International 
Symposium on Computer-Assisted Cartography - AUTO-CARTO III, pp. 
776-791, 1979. * 


1-8. General Electric Company, Preliminary Design Requirements for 
C ensus/Urbanized Area Applications Systems Verification and 
Transfer - Final Report , prepared under NASA contract NAS5-23412, 
Greenbelt, Maryland, 1977. 

1-9. General Electric (Space Division), Synopsis of Census-Urban Area 
ASVT Baseline Processing , Interim Report to NASA/GSFC, 1978. 

1-10. General Electric (Space Division), Test Phase Studies for the / 
Census-Urbanized Area Application System Verification and 
Transfer , Contract NAS5-23412, Mod. No. 106, 1978. !.! 


6-1 


Bryant, N. A., and A. L. Zobrlst, ’’Integration of Socioeconomic 
Data and Remotely Sensed Imagery for Land Use Applications, " 
Proceedings, Second Annual William T. Pecora Symposium , American 
Society of Photogramme try, pp. 120-130, 1976, 

Bryarit, N. A., and A, L. Zobrlst, "IBIS: A Geographic Information 

System Based on Digital Image Processing and Image Raster Data 
Type," IEEE Transactions on Geos cience Electronics , Vol. GE-X5, 

No. 3, pp. 152-159, 1977. 

Bryant, N. A., and A. L. Zobrlst, "An Image Based Information 
System to Correlate Imagery With the Users' Data Requirements," 
paper presented at the 30th Annual Conference of the Society of 
Photographic Scientists and Engineers, 1977. 

Friedman, S. Z., Use of the Ratio Threshold Classifier In Mapping 
Urban/Non-Urban Land Cover from Landsat MSS Imagery , Report No. 
715-25, Jet Propulsion Laboratory, Pasadena, Calif., February 1980 
(JPL internal document). 

NASA, Data Users Handbook: NASA Earth Resources Technology 

Satellite , Document 71SD4249, NASA Goddard Space Flight Center, 
Greenbelt, Maryland, 1972. 

Thomas, V. I., 1975, Generation and Physical Characteristics of 
the Landsat 1 and 2 MSS computer Compatible Tapes , NASA Goddard 
Space Flight Center, publication X-563-/5-223, 1975. 


Colvocoresses, A. P., Space Oblique Mercator: A New Map Projection 

of the Earth , Memo for the record August 6, 1974, by the Cartographic 
Coordinator, U.S. Geological Survey, Reston, Virginia, 1974. 

Kriegler, F. J. , et al., 1969, "Preprocessing Transformations and 
Their Effects on Multispectral Recognition," Proceedings, Sixth 
International Symposium on Remote Sensing of Environment , pp . 

97-131, 1969. 


Ready, P. J. and P. A. Wintz, "Information Extraction, SNR Improve- 
ment, and Data Compression in Multispectral Imagery," IEEE Trans- 
actions on Communications , pp. 1123-1130, 1973. 


Todd, W. J. , "Urban and Regional Land Use Change Detected by Using 
Landsat Data," Journal of Research , U.S. Geological Survey, pp. 
529-534, 1977. 


6-2 


0 


ft 

BIBLIOGRAPHY 

Anderson, R. , "Land :‘]se and Land Cover Changes - A Framework for Moni- 
toring," Journal of Research , U.S. Geological Survey, pp. 143-153, 


Angelici, G. L.„ and N, A. Bryant, "Techniques i for the Creation of Land 

Use Maps and Tabulations from Landsat Imagery," Proceedings , Second 
Annual William T. Pecora Symposium , American Society of Photo- 
grammetryj, pp. 162-171, 1976. 

Angelici, G. L. , and N. A. Bryant, "Solar Potential Inventory and 

Modeling, " Data Resources and Requirements: Federal and;Local 

Perspectives , Sixteenth Annual Conference of the Urban and 
Regional Information Systems Association, Washington, D.C.,pp. 
442-453, 1978. 

Borden, F. Y., and D. L. Williams. A Reduction in Ag. /Residential 
Signature Conflict Using Principal Components Data , GSFC 
Publication X92^/-77-182, Goddard Space Flight Center, 1977. 

Bryant, N. A., and A. L. Zobrist, "Integration of Socioeconomic Data 

and Remotely Sensed Imagery for Land Use Applications," Proceed- 
ings, Second Annual William T. Pecora Symposium , American 
Society of Fhotogrammetpy, pp. 120-130, 1976. ' 

Bryant, N. A., and A- L, Zobrist, 1977, "IBIS: A Geographic Information 

System Based on Digital Image Processing and Image Raster Data 
Type," IEEE Transactions on Geoscience Electronics , Vol. GE-15, 

No. 3, pp. 152-159, 1977. 

Bryant, N. A., and A. L. Zobrist, "An Image Based Information System to 
Correlate Imagery With the Users' Data Requirements," paper 
presented at the 30th Annual Conference of the Society of Photo- <- 
graphic Scientists and Engineers, 1977. 

Christenson, J. W. , J. B. Davis, V. J. Gregg, H. M. Lachowski, and 

R. L, McKinney, Landsat Urban Area Delineation , Intralab Project 
75-3, NASA/Goddard Space Flight Center, Greenbelt, Maryland, 1977. 

* ■; 4 / ° 

Christenson, J, W., et al, , -"Urbanized Area Analysis Using Landsat Data," 
paper presented at the Conference on Computer Mapping. Harvard 
University, Cambridge, Mass., 1978. O 

. : •. . ‘ 

Colvocoresses, A. P., Space Oblique Mercator: A New Map Projection of 

the Earth , Memo of the record August 6, 1974, by the Cartographic 
Coordinator, U.S. Geological Survey, Reston, Virginia, 1974. 


7-1 


Davis, J. B. and S. Z. Friedman, "Assessing Urbanized Area Expansion 
Through the Integration of Landsat and Conventional Data," 
Proceedings, 45th Meeting of the American Society of Photogrammetry , 
pp. 776-791, 1979. 

cf? 

Friedman, 3. Z., 1978, Use of the Ratio Threshold Classifier in Mapping 
Urban /Non-Ur ban Land Cover from Landsat MSS Ittagery , Report No. 
715-25, Jet Propulsion Laboratory, Pasadena, Calif., February 
1980 (JPL internal document) . < 

Friedman, S. Z., and G. L. Angelici, "The Detection of Urban Expansion 
from Landsat Imagery," Remote Sensing Quarterly , Vol. 1, J No. 1, 

1979; presented at Annual Spring Meeting of the Association of 
American Geographers, 1978. 

General Electric Company, Preliminary Design Requirements of Census /Ur- 
banized Area Applications Systems Verification and Transfer-Final 
Report , prepared under NASA contract NAS5-23412, Greenbelt, 

Maryland, 1977. 


General Electric (Space Division) , Synopsis of Census-Urban Area ASVT 
B aseline Processing , Interim Report to NASA/GSFC, 1978. 

General Electric (Space Division) , Test Phase Studies for the Census- 
Uibanized Area Application System Verification and Transfer , 
Contract NAS5-23412, Mod. No. 106, 1978. 

Kriegler, F. J. , et al. , "Preprocessing Transformations and Their Effects 
on Multispectral Recognition," Proceedings , Sixth International 
Symposium on Remote Sensing of Environment ," pp. 97-131, 1969. 

McKinney, R. L. and Stauffer, M. L., Landsat Image Differencing as. an 
Automated Land Cover Change Detection Technique , Computer 
Sciences Corporation/TM-78/6215, 1978... 

NASA, Data Users Handbook: NASA Earth Resources Technology Satellite , 

Document 71SD4249, NASA Goddard Space Flight Center, Greenbp.lt , 
Maryland, 1972. 

NASA, Project Plan: Applications Systems Verification Test (ASVT) , 

Census Urbanized Area , NASA Goddard Space Flight Center, Greenbelt, 
Maryland, 1976. o ? 


7-2 


Ready, P. J, , and P, A. Wintz, 1973, "information Extraction, SNR 
Improvement, and Data Compression in Multispectral Imagery," 

IEEE Transactions on Communications, pp. 1123-1130, 1973. 

Seidman, J. B., and A. Y. Smith, VICAR Image Processing System: Guide 

to System Use , JPL Publication 77-37. Revision 1, Jet Propulsion 
Laboratory, Pasadena, Calif., 1979. 

Thomas, V. I., Generation and Physical Characteristics of the Landsat 
1 and 2 MSS computer Compatible Tapes, NASA Goddard Space Flight 
Center , publication X-563-75-223, .1975 . 

Todd, W. J., "Urban and Regional Land Use Change Dete/ ' . d by Using Land- 
sat Data," Journal of Research , U.S. Geological * i'vey, pp. 529- 
534, 1977. ’ V 

o 

Wherry, D. B., and S. Z. Friedman, 1979, "Cartographic Applications of 
an Image Based Information System," Proceedings , International 
Symposium on Computer-Assisted Cartography - AUT0-CART0 III, pp. 
776-791, 1979. 

Yap, F. I., The Application of an Image Based Information System (IBIS) 
to Urban Classification Studies — Test Case: Richmond, Virginia , 

Computer Sciences Corporation/TM-78/6196, 1978. 

!i 



APPENDIX A 


A TECHNICAL DESCRIPTION OF THE 


IMAGE-BASED INFORMATION SYSTEM 

In Section II the Image-Based Information System (IBIS) was 
described on a conceptual level. A more concise and technical 
description of IBIS is included in this appendix. Most computer „ 
programs comprising the system are described in summar y J 30 ?£ Annua i 
following article is a facsimile of a paper presented at fche 30th Annual 
Conference of' the Society of Photographic Scientists and Engineers, 

Los Angeles, May 1977. 
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AN IMAGE BASED INFORMATION SYSTEM: 
ARCHITECTURE FOR CORRELATING SATELLITE AND TOPOLOGICAL 

DATA BASES 


I 

N.A. Bryant and A.L. Zobrist 


INTRODUCTION 

Geographic information systems should satisfy four basic criteria if 
tuoy are to be useful: 1) They should provide specific point locations, 

as well as area locations of data; 2) They should provide for variable 
aggregation (sub-setting) of data; 3) They should provide a method for 
representing spatial arrangements; and 4) They should be able to interface 
with mathematical and statistical programs which can be called as needed 
to aid in the analysis of spatially-oriented data. * Practitioners of the 
art of geocoded systems design have progressed with varying degrees of 
success towards the goals outlined. As a rule, generalized systems have 
only rudimentary data manipulation capability (i.e., status updating and 
interrogation by area) , while highly specific and specialized systems have 
progressed further with modelling applications." 

In response to the desire to access data for selected areas, polygon 
and grid cell geocoded information systems have been developed. Such 
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systems rely on the tabular formatting of the input data, a costly and 
tine consuming process. Often the system falls into disuse because the 
updating of major segments of the data base becomes prohibitively expensive. 
In response to the desire to provide up-to-date resource infoimation, 
investigators have studied the feasibility of applying Landsat and high 
altitude photo imagery to natural resource mapping. The consensus which 
appears to be evolving is that , while remotely sensed imagery can provide 
timely coverage and sufficiently accurate maps using ADP techniques, the 
end product is still a map that cannot interface directly with an existing 
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geocoded information storage and retrieval system. 

FUNCTIONAL REQUIREMENTS « 

Data Management Considerations 

The ease with which an agency can establish a geographic information 
system is constrained by the level of detail and the computer technology 

available. With geocoded data, there is a dramatic increase in file size 

■ ., ■ • »> 

with every added variable and increased resolution. Thus, a parcel level 
system for the City of Los Angeles consists of over 800,000 records, while 
a census tract level system has less than 1,000. The use of newer genera- 

•• .. o 

tion computers has only moderately improved the overall operation, as the 
major efforts are involved in both file generation and editing and computer 
software architecture. Frequently, urban information systems have become 
underutilized because of high' manpower requirements for the initial encoding 
of data and the heavy expense involved in updating. Land use character- 
istics have consistently been the most expensive data to incorporate in 
information systems, primarily because there has been no .subsidy available ? 


to the user in recording and preparing resource inventory data, such as 

4 

occurs with demographic, economic, health, and assessor data sets. 

There are a number of ways to encode spatial data . 5 It is possible, 
however, for illustration purposes, to dichotomize referencing systems into 
nominal and ordinal, and data types into tabular, graphical, and image. 
Ordinal systems reference data by the actual geographical coordinate values. 
Thus, natural resources such as forests, rivers, crop lands, and geologic 
formations are mapped with selected identifiers (total area, boundary, 
centroid) referenced to latitude and longitude or other selected geographic 
coordinate system. Nominal systems are "name referencing", i.e. , data or 
information is referenced to a name-designating system. Any district- 
based referencing convention, such as census tract, sewer district, town- 
ship, or transportation zone, assumes the operator knows where each 
administrative area is located and leaves the analysis of contiguity- 
effects, etc., up to the individual.^ Problems invariable arise during 
the computational processing involved in conversions between ordinally and 
nominally referenced data and have, with the exception of the Census DIME 

(Dual Independent Map Encoding) • File system, forced the conversion of all 
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data to one reference system or the other. 

Of equal concern to the designer of a geographic information system 
is an incorporation of the various geocoded data types to assure adequate 
"data capture". It is only through the integration of tabular, graphical, 
and image data types that geographic information systems achieve the 
syneTgistic impact required to offset their initial cost. For instance, 
tabular tabular files may keep records of individual weather stations, 
while graphical files record elevation contours, and image data sets the 


distribution of land use. The combination of all three data types would 
provide analysts with the variety of spatial data needed to model alterna- 
tive sources of non-point source air pollution. The need to deliver a 
uniformly encoded result to the user has been a key element to changes 

that have occurred in both geocoding approaches and computer system archi- 

o 

tecture applied to geographic information systems. 


Approaches to Geocoding 

There exist three principal goecoding architectures, which have evolved 
from simple grid cell systems to the complex polygon methodologies, and 
most recently include the image raster data type (see Figure 1). 


METHOD 



COST CONSIDERATION 
MANUALLY OPERATED 


OVERALL FLEXIBILITY 


SPATIAL RESOLUTION 
POOR, UPDATING 
DIFFICULT 


POLYGON: 



EXPENSIVE FOR URGE 
DATA SETS 


REQUIRES IMAGE 
PROCESSING TECHNOLOGY 


CERTAIN OPERATIONS 
PROHIBITED 


NEITHER SCALE NOR 
DATA FORMAT DEPENDENT 


FIGURE 1. Approaches to Geocoding 

From several recent reviews, it is evident that the goals of geographic 
information systems have become more ambitious as faster computer systems 
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evolved and peripherals become more sophisticated. Grid cell methods 
serve the need to retrieve geo- located data and generate maps through the 
cross tabulation of variables encoded within a particular cell. Several 
important drawbacks reduce the overall flexibility of grid cell systems: 
a) Their spatial resolution is only as accurate as the grid cell site 
(usually ranging from one acre to ten square miles) ; b) The systems permit 
the referencing of data in either a nominal or ordinal manner, never both; 
c) The need for manual encoding of the input data files has made updating 
difficult and- even prohibitively expensive and effectively limited the 
spatial resolution of grid cells to satisfy the need to achieve regional 
coverage. 

In response to the failings of grid cell geocoding, polygon systems 
grew as electronic coordinate digitizers became generally available. Poly- 
gon geocoding formats effectively solved the spatial resolution dilemma 
inherent with grid cell formats, while coordinate digitizing hardware has 
permitted rapid encoding of data. The most important achievement, however, 
was the integration of nominal and ordinal referencing in the Census DIME 
file methodology. 1 ^ Despite these significant achievements, polygon geo- 
coding systems have left the problem of ordinal data updating unresolved 
and created new challenges inherent in their graphical data structure. 

These problems include: a) considerable conputational expense associated 

with file editing; b) complex topological architectures to achieve efficient 
data extraction from any given area;" 1 and c) large investments ' in computer 
systems to achieve polygon overlay of separate files for encoding ordinal 
data into nominal encoding formats (e.g., acreage of land use for each 
census tract) . Many of these constraints can be mitigated, by the use of 
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an image raster encoding procedure and application of digital image process* 
ing algorithms to implement geographic information system analyses. 

IBIS Data Management 

o 

Digital image processing techniques can be applied to interface 
existing geocoded data sets and information management systems with thematic 
maps and remotely sensed imagery. The basic premise is that geocoded data 
sets can be referenced to a raster scan that is equivalent to an ultra* 
fine mesh grid cell data set, and that images taken of thematic maps or 
from remote sensing platforms can be converted to a raster scan. Until 
recently, the image format has been used primarily as a computer proces sable 
equivalent of a photograph, with the value stored in each cell of the image 
representing a shade of grey or a color. But if the image is of a geo* 
graphical point can be accessed immediately by position in the image matrix. 
Figure 2 illustrates the calculation of memory address of the data value 
from a latitude- longitude pair. 

The image datatype seems to be a powerful and general representation 

for spatially distributed data, and the range of uses can be divided into 

l 7 

several broad categories: 

Physical Analog : The pixel value represents a physical variable such as 
elevation, rainfall, smog density, etc. 

District Identification : The pixel value is a numerical identifier for q 
the district which includes that pixel area. 

Class Identification : The pixel value is a numerical identifier for the 
land use or land cover, or for any other area classification scheme. 
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LOCATION FORMULAS: 

P 

*.[*«■ UI *1 . LONG ♦ e ] nearest integer 
Y- [d.LAT.E » LONG + f] NM , Er , Nt(sf|t 

MEMORY ADDRESS • BASE + KX + Y 

FIGURE 2 Image Matrix As A Data Representation 
Tabular Pointer : The pixel value is a record pointer to a tabular record’ 

which applies to the pixel geographical area. 

Point Identification : The pixel value identifies a point, or the nearest 

of a set of points, or the distance to the nearest set of points. 

Line Identification : The pixel value identifies a line, or the nearest of 

a set of lines, or the distance to the nearest of a set of lines. 

This range of uses required that the system be able to handle images 

composed of words of varying length. For example, to identify census tracts 

« » 

in Los Angeles County requires 1500 different pixel values and elevation 
maps can require 15,000. This is more than the usual 256 grey levels 
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handled by photographic image processing systems. 


For each data type, the overriding consideration is usually getting 
the data into the system for it is in this area where greatest costs and 

o 

difficulties are usually incurred. Figure 3 depicts data input as a three 
stage process. The first stage, called data capture, includes all opera- 
tions up to the point where a data file is computer readable. Data capture 
costs are enormous for many basic kinds of data, for example, demographic 
and economic data gathered by the U.S. Bureau of the Census. These data 
are then made available on computer tape at nominal cost to any user. 

•') V 

Another common method of data capture is coordinate digitization of 
boundaries or linear features from a map. The map is not computer compatible 
but the digitizer output is. Manual photointerpretation (for example, 
manual determination of land use boundaries) is a common step prior to 
coordinate digitizing. Landsat imagery is a particularly attractive data 
source because it is already computer conpatible and gives up-to-date 
coverage of large areas at nominal cost. Editing costs vary widely 
depending upon the nature of the data source. Reformatting can also be 
a major operation where large data files are involved. A good example is 
community analysis where data are gathered by various districtings (police, 
fire, sewer, census tract, etc.) but must be reformatted to a coninon 
districting so that analysis can be performed (for example, to obtain 
police calls per capita). The final stage of data input is to obtain a 
temporal baseline. For example, police calls in 1976 cannot be divided 
by 1970 population to obtain a per capita figure because of a temporal 
difference. Establishing a temporal baseline involves projection and 
modelling of the initial data to one or more coninon points of time. 
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FIGURE 3 Data Input Stages 

IBIS Data Analysis 

Once a working data base is set up, provisions must be made for 
information retrieval, information analysis, and report generation. 
Operations here are usually of a much smaller scale than data capture in 
terns of time and cost, but there is instead a question of flexibility 
and ease of use, and a question of system cost. All data analyses can be 
laid out as a sequence of primitive steps, thus, a functional requirement 
is that an adequate set of primitive operations can be implemented. 
Mathematical and statistical analyses of tabular files are well -understood, 
and packaged systems can easily be interfaced. The open question is 
whether geo-based file computational steps can be implemented. Some 
examples of these are: 1) Given a point and a district, does the point 

lie within the district. 2) Given a point and a district file, which 



district contains the point. 3) Given a particular district in a district 
file, what are its neighbors. 4) Given a district file and an area 
classification file, what are the acreages of each area classification in 
each district. 5) Given two district files, one major and one minor, what 
are the proportions of each minor district in each major district, 6) 

Given a district file and a line segment, what are the mileages of the 
line segment in each district. 7) Given a point p and a point file, what 
is the point in the point file nearest to p. 8) Given a point p and a 
point file, what is the distance from p to the nearest point in the point 
file. 9) Given a point and a line segment file, which line segment passes 
closest to the point. 10) Given a density map and a district file, what 
are the volumes in each district (spatial integral of density). 11) Given 
a district, what is the centeimost point (an inside point which is farthest 
from the boundary). And for the connoisseur: 12) Given a district file, 

assign four color? to the districts so that a map can be produced with 
adjacent regions always a different color. 

.The preceding list is just a sample of the sorts of spatial or geometric 
calculations which need to be performed by a comprehensive geo -base infor- 
mation system. More complex operations will usually be implemented as a 

sequence of these primitive operations, but because of the magnitude of 

jj \ ■ 

the data files' and because of iteration due to modelling, compute time can 
// ■ 
be a serious problem. A method which solves one of the primitive problems 

in 0.1 second may seem usable, but not if it has to be performed ten million 

times for a particular application. 

It is worth noting here that many of these operations are difficult 
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and time consuming if the 'working data base is in polygon or graphical 
format (i.e.# lines &re specified by their end points and a district is 
given by a sequence of line segments). In particular, the operation called 
polygon overlay which solves primitive problems four and five is extremely 
difficult to perform on large files in graphical format. If the files are 

in image format, then polygon overlay becomes a simple counting operation 

h • 17 \ . 

(Figure 4). > ^ 

/! , ‘ 

A different area of concern is the interfacing of nominal and ordinal 
geo-referenced files. Given that part of the working data base has been 
obtained (or reformatted) into an ordinal format, the major part of the 
data base will probably still be in tabular fori..* As an example, an air 

pollution density image might be interfaced with a population table to 

>1 ■ 

obtain a measure of health effect. Speaking more generally, the system 
must be able to perform operations on mixed data types, thus allowing the 
working data base to be built from raw data in its most natural form. 

An operation of special importance is crosstabulation, mentioned 
previously in connection with data management. It converts data aggregated 
by one district convention to an aggregation by another district convention. 
The operation itself is trivial since it involves multiplying by a set of 
factors which measure the percentage subareas of one districting in the 
other. The crosstabulation factors are derived by polygon overlay and may 
also be modified by a density estimate of the variable being crosstabulated. 

O) ' . . . : •" _ _ ' ' '■ ■ ' ' 

A comprehensive system should be able to represent district data sets in 
such a way that they can easily be updated and the factors rederived. This 
means that the tabular data can be kept in its natural unit of aggregation 
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in the working data base and quickly crosstabulated to whatever districting 


is needed as analysis proceeds. 
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GRAPHICAL FORMAT: 

0 


DISTRICT 


DISTRICT 

FILE 1 


FILE 2 

(M DISTRICTS) 


(N DISTRICTS) 


MUST CALCULATE 
AREA OF M x N 
INTERSECTIONS 
OF POLYGONS 


IMAGE FORMAT: 



■ <;;> 

MUST COUNT 
K x L PIXEL 
PAIRS 




FIGURE 4 Comparison of Polygon Overlay Techniques 


Because the image datatype 1 is used, capabilities for digital image 
file handling, image manipulation, and image processing are required. Thus, 

. . ■ . ' O'"-’: . 

the IBIS system has been built upon an existing image processing system, 

13 

VICAR (Video Image Communication and Retrieval), developed at JPL. 

Certain basic image processing operations are absolutely essential. One 
must accomplish image- to- image registration, whereby images of different 
scale, rotation, or map projection are superimposed precisely enough so 
that corresponding pixels represent the same geographic location. Rubber- 
sheet registration is almost always necessary to achieve, the needed degree 



it 
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of accuracy. On other hand, it is anticipated that even esoteric image 
processing operations, such as convolution smoothing will be useful for 
certain types of applications. The conclusion here is that any image - 
based information system must contain a powerful image processing subsystem. 

Finally, the system must be capable of computations which will allow * 

. '■ : x ■ ■ I 

it to tun a variety of models. The key here is the introduction of the 



image raster as a spatial data representation. Contiguity effects can 


then be handled because adjacent areas are contiguous in the image raster 
(in terms of their computer address) whereas, with other file representa- 
tions, adjacent areas are accessed by file searching. Basic operations. 


such as spatial integration with reporting by district are greatly simplified \ 

... , - , ' ■ ■ $ ■ / T - 1 

with an image representation. More esoteric models, say, involving cellular j 

o % I 

transformations or diffusion processes may already use a matrix on raster j 

\ representation but will be helped by the synergism of placing their j 

specialized capabilities into a general and comprehensive data management \ 

system. 0 q 

.. . r ' O . 

SYSTEM DESCRIPTION 

IBIS is presently resident on an IBM 360/65 system with one megabyte 
of core with interfaces to a large number of disks and tape drives to 
accommodate the large files associated with image processing technology. 

The operating system is OS 360/MVT and all IBIS routines are designed to 
run in a 150 kilobyte region. Present applications use less than one 




twentieth of the total machine time, so this, hardware configuration would 
be suitable for a large scale data processing operation at the national 
level. A mini-computer version is partially complete. Although the total 

A-14 


S"-‘ 


cost of an IBIS system Is high, the unit cost of operation is low, 
especially when compared to existing methods of performing the types of 

tasks the system is intended for. 

the user request is given to the IBIS system by means of a language 
which is translated into the host machine job control language. The trans 
lated code can then invoke system functions or processing modules. This 
organization makes the system flexible and easily extendable. The system 
software consists of a number of FORTRAN modules and a relatively small 
nucleus. All of the factors mentioned here (modular design, use of 
FORTRAN, good interfaces to user and hardware) made technology transfer 
more feasible. 

The processing modules constitute a primitive set of functions 
operating on the various datatypes to achieve the functional requirements 
set out in the previous section. Figure 5 is a schematic layout of the 
central part of the IBIS system, emphasizing the interfaces between image, 
tabular, and graphical datatypes. 

The first set of routines convert graphical (polygon) data files to 
digital image descriptions of regions or areas. 

VTRACTI, SGRBGEN, VSCR3GEN - All three programs basically have the same 
function, convert polygon files described by x,y coordinates to a VICAR 
record format, but each uses a. different input type: 

VTRACTI - Census Tracts 

SCRBGEN - Cards 

VSCR3GEN - Tppe 
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GRAPHICAL 

DATA 

FILE 
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POLYREG 

POLYGEOM 

POLYSCRIBE 

FILL 

PAINT 


IMAGE 

PROCESSING 

(VICAR) 


( POLYOVLY 
( TALLY 

( SORT 
AGGRG 
AGGRG2 
MF 

CROSSTAB 

TRANSFER 

ZIPCOL 

COPIN 

COPOUT 


FIGURE 5 IBIS Schematic 

POLYREG - Performs rigid rotation and scaling of polygons prior to input 
into POLYSCRB.) 

POLYGEOM - Performs rubber sheet mapping on polygons prior to input into 
POLYSCRB. r- ■ 4 , k . . ■ '/' - ' i ■ 

POLYSCRB - Scribes lines into an image, file, by setting the boundaries of 
polygons to a particular DN (grey scale) value. Tne routine was 
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especially designed for the conversion of thousands of lines and has 
parameters for chaining lines and closure of polygons. 

FILL - Fills holes, thickens lines and/or removes noise with variable 

thresholds and window sized for enhancement of scribed polygons prior 
to PAINT. 

PAINT - Converts an image with scribed polygons into a multi-color map where 
the DN corresponds to a map color. The routine will handle up to 
30,000 polygons or regions. 

Moving from top to bottom in Figure 5, the VTCAR image processing 
routines are next. The VICAR system presently has over 100 routines for 

the manipulation of images of varying Work length. Of special interest to 

I 

modellers is the following routine. 

f 

F2 - Performs array arithmetic on a pair of images in byte or halfword 
format. The function to be applied is specified by a FORTRAN-like 
expression. 

Polygon information extraction routines perform the overlay operation 
on two sets of polygons, the reference polygon created in paint and digital 
image data. Each pixel on the image file is matched with the corresponding 
pixel on the reference file and stored by pixel pairs. The following 
programs are used in the overlay process. 

PQLYOVLY - Produce a histogram (or pixel count) or registered images by 
DN values for single pixels or a joint histogram by DN pairs as 
generated in the overlay procedure. Since the problem of storage- is 
critical to this process, the user has three optional storage methods 
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specified in the parameters. 

TALLY- Produces a histogram similar to POLYOVLY except that the DN 
values in one of the two images are sunned. 

The results of POLYOVLY and TALLY generate data files stored by 
column. The colunn organized file is called the interface file. This 
file serves as an interface between image and tabular data sets. The 
interface file is manipulated to produce a report by the following routines: 

SCRT - Sorts the interface file into ascending or descending order 

# 1 . 

according to one or more columns. 

AGGRG - Aggregates columns of numbers using a designated column as an 
index for summation in other columns. 

AGGRG2 - Summarizes and collates columns of numbers to produce a single 
row for each control value in the column which is used as an index. 

MF Performs column arithmetic given an arithmetic expression in the 
parameter field which denotes columns such as: 

"CS-C2/C3+100+SQRT (C2) " . 

CROSSTAB - Tabulates information referenced to one polygon districting to 
another polygon districting through the polygon overlay technique. 
TRANSFER - Change vertically aligned columns of data as produced in 
POLYOVLY and used in previous routines to smaller vertical columns 
based upon data values (e.g., land use). 

ZIPCOL - Substitutes column index values with user district names or 
numbers. . -W 

After processing by the above routines, the interface file is ready 
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for report generation. 


MPRINT - Prints out the interface file according to a relatively simple 
format. 

REPORT - Prints out the interface file with user specified titles, 

column titles, paging, spacing, etc. Subtotaling on a given control 
column can be requested. There are numerous features for formatting 
and for printing alphabetic data. 

To interface tabular files and also to give tape output as an 
alternative to report generation,, the following routines are used. 

COPIN - Copies columns of a tape file into columns of the interface file 
according to an index column in the interface file. The index must 
be present in the tape file. The files must be sorted according to 
the index and a merge is performed. 

COPOUT - The reverse of COPIN. Columns are written from the interface 
file to the tape-. 

CONCLUSION; 

An image-based information system is necessary for the full utilization 
of satellite imagery data and anticipated development of regional land 
capability analysis centers. The future availability of frequent updates 
of land resource inventory statistics, with a known and acceptable sampling 
accuracy, should permit the incorporation of this data with the annual up- 
dates published by other governmental bureaus. An image based information 
system does more than introduce remotely sensed imagery t.o the mainstream 
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of data processing application; it provides a new approach to the management 
and analysis of spatially-referenced data. 

The projected demands tc; be placed upon geographic infonnation 
systems will place a strong emphasis upon the capability to store and 
retrieve large amounts of data and manipulate data sets for portions of 
the files efficiently. A major drawback facing most geocoding procedures 
is that they rely on sequential computations applied to tabular data 
strings, and as such require a large investment in formatting or process- 
ing data that is inherently two-dimensional. Raster scan data bases will 
avoid many of these problems and possess additional advantages. The video 
communications field has been and is continuing to address both the problems 
of mass storage and applications of rapid interactive processing that 
places a minimal reliance upon computer software routines. The specialized 
requirements of geographic information systems should derive considerable 
benefit from the image processing field in the future. 
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APPENDIX B 

SELECTED STATISTICAL REPORTS 

The tabular reporjjt summarizing Information stored in an IBIS 
data base is a major feature of the Image-Based Information System, 
These reports may be useful to the Census Bureau geographer in de- 
termining areas of urban expansion. 

Two reports are included in this appendix. The first report 
is a complete summary of urban land cover in 1970 and 1975 for the 
Orlando SMSA. The second report is a limiteid summary of land cover 
in 1970 and 1975 for the Seattle-Everett SMSA. Only those tracts 
that were found to be nonurban in 1970 have been listed. 



CENSUS 

TRACT 

number 


101.00 

102.00 

103.00 

104.00 

105.00 

106.00 

107.00 

lot). 00 

109.00 

110.00 
111.00 
112.00 
111.00 

114.00 

115.00 

116.00 

117.00 

118.00 

119.00 

120.00 
121.00 
122.00 

123.00 

124.00 

125.00 

126.00 

127.00 

128.00 

129.00 

130.00 

131.00 

132.00 

133.00 

134.00 

135.00 

136.00 

157.00 

138.00 

139.00 




Table B-l. 

URBANIZED LAND COVER STATISTICS FOR THE ORLANDO SNSA 
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1970 AND 1975 
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40 

140.00 
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0 
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41 
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42 
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0 
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100.0 

YES 
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0 

0.0 

146 

2.9 

146 

100.0 


5l 

151.00 

5075 

6331 

798.4 

0 

0.0 

130 

2.6 

130 

100.0 


52 

152.00 

2762 

64 79 

1501.5 

2762 

100.0 

2 762 

100.0 

0 

0.0 


53 

153.00 

950 

4285 

2887.5 

950 

100.0 

950 

100.0 

0 

0.0 


°54 

154. 00 

1187 

4936 

2661.6 

1187 

100.0 

1187 

100.0 

0 

0.0 


55 

155.00 

1897 

3468 

1169.9 

1897 

100.0 

1897 

100.0 

0 

0.0 


56 

156.00 

1000 

4908 

3142.5 

1000 

100.0 

1000 

100.0 

0 

0.0 


57 

16/ • Oil 

952 

3181 

2138.2 

952 

100.0 

952 

100.0 

0 

0.0 


68 

158.00 

1032 

4286 

2659.1 

1032 

100.0 

1032 

100.0 

0 

0.0 


59 

169.00 

717 

4361 

3590.7 

717 

100.0 

717 

100.0 

0 

0.0 


60 

160.00 

964 

4216 

2799.0 

964 

100.0 

964 

100.0 

0 

0.0 


61 

161.00 

1142 

4355 

2441.0 

1142 

100.0 ' 

-1142 

100.0 

0 

0.0 


62 

162.00 

804 

6235 

4963.7 

804 

100. 0 

804 

100.0 

0 

0.0 


63 

163.00 

1157 

6162 

3407.8 

1157 

100.0 

1157 

100.0 

0 

0.0 


64 

164.00 

7420 

3959 

341.5 

0 

0.0 

687 

9.3 

687 

100.0 


61> 

lu5.00 

9319 

6067 

279.3 

0 

0.0 

0 

0.0 

0 

0.0 


66 

166.00 

88141 

2224 

16.1 

0 

0.0 

6 

0.0 

6 

100.0 


67 

167.00 

157858 

7379 

29.9 

0 

0.0 

1483 

0.9 

1483 

100.0 


68 

168.00 

41116 

6369 

99.1 

0 

0.0 

2864 

7.0 

2864 

100.0 


69 

169.00 

3943 

5371 

871.9 

0 

0.0 

1374 

34.9 

1374 

100.0 

YES 

70 

170.00 

24780 

3345 

86.4 

0 

0.0 

1351 

5.5 

1351 

130.0 


71 

171.00 

76679 

3817 

31.9 

0 

0.0 

2406 

3.1 

2406 

100.0 


72 

172.00 

1749 

1007 

368.5 

0 

0.0 

0 

0.0 

0 

0.0 


73 

173.00 

1240 

2911 

1502.1 

1240 

100.0 

1240 

100.0 

0 

0.0 


74 

174.00 

2564 

4811 

1205.6 

2554 

100.0 

2564 

100.0 

0 

0.0 


76 

175.00 

14725 

5823 

263.1 

0 

0.0 

£327 

9.0 

1327 

100.0 


76 

176.00 

787 

2633 

2140.3 

797 

100.0 

787 

100.0 

9 

0.0 


77 

177.00 

2795 

3197 

732.1 

0 

U.O 

72 

2.6 

72 

100.0 


78 

178.00 

30064 

4600 

98.0 

0 

0.0 

1694 

5*6 

1694 

100. 0 
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Table B-l. (continuation 2) 


CENSUS 

TRACT 

NUMBER 


TOTAL 

ACRES 

PER 

TRACT 


1970 POPULATION 
STATISTICS 


NUMBERS DENSITY 
OF PER MILE 
PEOPLE SQUARED 


4970 URBANIZED 
LAND COVER 
STATISTICS 


ACRES PCT 


1975 URBANIZED 
i«NU COVER 
STATISTICS 


ACRES 


““JSSit UReIS*" 

4970 ANO 1975 


ACRES 


PCT MAJOR 


79 

179.00 

30379 

80 

201.00 

665 

81 

202.00 

596 

82 

203.00 

732 

83 

206.00 

1359 

86 

205.00 

667 

85 

206.00 

6369 

86 

207.00 

21553 

at 

208.00 

12696 

88 

209.00 

3130 

89 

210.00 

11803 

90 

211.00 

3565 

91 

212.00 

66936 

92 

213.00 

26396 

93 

216.00 

6035 

96 

215.00 

5227 

95 

216.00 

11830 

96 

217.00 

6705 

97 

218.00 

2356 

98 

219.00 

2077 

99 

220.00 

2633 

100 

221.00 

6526 

lOl 

222.00 

5517 

102 

NO TRACT 

1161155 


1960873 


3533 
2862 
3266 
3253 
6750 
3223 
1661 
1637 
3889 
5676 
963 
2 865 
2131 
3620 
1650 
6960 
3519 
5620 
7036 
3766 
7913 
3865 
5693 
0 


628009 


76.6 

6090.6 

3681.1 

2865.8 

2237.2 

3186.2 

215.0 

68.6 

199.2 

1201.6 

51.1 

513.6 
21.0 
87.8 

175.0 
60*9 

190.6 

766.5 

1910.8 
1153.5 
1923.1 

566.5 

660.6 

0.0 


0 

665 

596 

732 

1359 

667 

0 

0 

O 

3130 

0 

0 

0 

0 

0 

0 

0 

0 

2356 

2077 

2633 

0 

0 

0 


0.0 

100.0 

100.0 

100.0 

100.0 

100.0 

0.0 

0.0 

0.0 

100.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

100.0 

100.0 

0.0 

0.0 

0.0 


77365 


2026 

665 

596 

732 

1359 

667 

296 

671 

1383 

3130 

686 

0 

0 

0 

636 

285 

1253 

689 

2356 

2077 

2633 

726 

680 

0 


105966 


6.7 

109.0 

100.0 

100.0 

100.0 
100.0 

6.8 

3.1 

il.l 

100.0 

5.8 

0.0 

0.0 

0.0 

l i:l 

10.6 

10.6 

100.0 

100.0 

100.0 

16.0 

8.7 

0.0 


100.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

100.0 

100.0 

0.0 

100.0 

0.0 

0.0 

0.0 

100.0 

loe.o 

lOv/.O 

100.0 

0.0 

0.0 

0.0 

100.0 

100.0 

0.0 
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APPENDIX C 

J ^ 

COMPUTER PROCESSING TIMINGS 

\Q 

Three tables containing timings of major computer processing 
steps for the three Census-Urbanized Area ASVT applications have 
been compiled. Timings refer to data processing on an IBM 360-65 
located at JPL. Both the elapsed time (wall-clock time) and the 
compute time (central processing unit time) have been tallied in 
minutes of execution. 

The table summarizing the Orlando application is most detailed, 

.0 

and has been divided into two segments, IBIS and VICAR processing 
steps. Processing summaries for the Seattle-Everett and Boston appli- 
cations have been listed in the sequence of execution. 


C-l 




*K==K#Wii 


" o ; . .?<Co - 

... f \ 

t 

Table C-l. COMPUTER PROCESSING TIMINGS ORLANDO SMSA APPLICATION 


IBIS Processing Steps 0 

General surface fit qf census tract file to planimetric base . . 

Final adjustments to complete fit (rubber-sheet!^.*') 

Region identification (formation of the geo-reference plane) ....... 

Process census tract centroid data ..- .......... 

Polygon overlay of geo-reference plane (to obtain 1970 area calculations) 

Generate 1970 population statistics 

Generate 1970 urbanized area map \ 

Generate 1970-1975 urban expansion map (data plane integration) 

Polygon overlay of expansion map and geo-reference plane (for statistics) 

Derive statistics and report (urban expansion 1970-1975) 

TOTAL IBIS PROCESSING TIME 


VICAR Processing Steps 


Logging of 1973 Landsat imagery (planimetric base) ... 

Logging of 1975 Landsat imagery (for 1975 urbanized land cover) . . 

Extraction of test area from 1973 imagery 

Registration of 1975 imagery to 1973 Imagery (rubber-sheeting) ... 
Eigenvector transformations of 1973 and 1975 imagery ........ 

Mask water bodies from eigenvector 2 of both 1973 and 1975 imagery . 
Mapping of urban/non-urban land cover for both 1973 and 1975 imagery 

TOTAL VICAR PROCESSING TIME .... 


TOTAL PROCESSING TIME 


TIMING IN MINUTES 
ELAPSED TIME COMPUTE TIME 

15.39- 

0 

2.87 

11.40 

4.48 

32.99 

11.22 

11.17 

2.05 

6.71 

1.17 

1.58 

0.25 

39.41 

16.55 

23.31 

3.13 

15.05 

2.38 

5.32 

1.05 

162.33 

(2:42.33) 

45.17 

(0:45.17) 

TIMING IN MINUTES 
ELAPSED TIME COMPUTE TIME 

131.28 

52,35 

79.10 

o ■ 56.50 

11.80 

1.97 

(y 

40.04 

7.98 

126.35 

30.57 , 

16.80 

4.62 

24.73 

4.75 

=^430.10 

(7:10.10) 

158 * 74 
u (2:38.74) 

592.43 

(9:52.43) 

203.89 

(3:23.89) 




Table C-2. COMPUTER PROCESSING TIMINGS SEATTLE-EVERETT SMSA APPLICATION 


0 

1 

w 


1. General surface fit of census tract file to planimetric base ...... 

2. Einal adjustments:-. to complete fit (rubber-sheeting) 

3. Region identification (formation of the geo-reference plane) ...... 

4. Process census tract centroid data 

5. Polygon overlay of geo-reference plane (to obtain 1970 area calculations) 

6. Generate 1970 population statistics * 

7. Generate 1970 urbanized area map . 

8. Extract 1975 land cover data from larger classification image . . . . . 

9. Generate 1970-1975 urban expansion map (data plane integration) c ... . 

10. Polygon overlay v Expansion map and geo-reference plane (for statistics) 

\ \ 

11. Derive statistics ana report urban expansion 1970-1975 ....... o» • 

TOTAL PROCESSING TIME 


TIMING IN MINUTES 

ELAPSED TIME 

COMPUTE TIP 

18.21 

3.88 

146.85 

56.77 

79, fa. 

31.28 

57.71 

29.64 

3.34 

2.12 

1.48 

0.31 

41.83 

12.56 

5.96 

3 ' 1-40 

62.00 

22.09 

\ 28.05 „ 

6.35 

5.58 

1.62 

450.60 

166.62 

(7:30.60) 

0 (2:46.62] 
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Table C-3. COMPUTER PROCESSING TIMINGS BOSTON REGION APPLICATION 
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NASA— JPL— Coml., L.A., Calif, 




